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Abstract 

We study the renormalization of the Wilson loop with a path corresponding to the Drell-Yan lepton 
pair production in two-loop approximation of perturbation theory. We establish the renormaliza¬ 
tion group equation in next-to-leading order and find a process specific anomalous dimension Toy 
in the corresponding approximation. 
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1. Introduction. The increase of the experimental accuracy necessitate the construction of the 
theory of power corrections for hardronic reactions. For the time being there is a firm theoretical 
ground only for the processes which admit the operator product expansion. However, these are 
the processes without the latter which are of the main interest. An important example of such 
a process is the Drcll-Yan (DY) lepton pair production with invariant mass Q 2 . Recent studies 
apparently reveal the correspondence between the soft-gluon resummation which is necessary for 
the processes going near the boundary of the phase space Jl], |J 0] and the non-perturbative power 
corrections 1,0- It was found that the ambiguity in the resummation of soft gluons to the leading 
logarithmic accuracy manifest itself in the power-like behaviour in the hard momentum scale lead¬ 
ing to power correction of the order Ai/Q ||. It was argued that this behaviour can be associated 
with particular matrix elements of Wilson line operators [0. The nature of this ambiguity was 
questioned in Ref. || where the authors attempted to identify the source of linear term in the 
perturbative series for the Wilson loop corresponding to the DY production. It was found that 
soft gluon resummation which allows to apply the renormalization group methods and infrared 
renormalons are disconnected since the former is related to the convergent analytical anomalous 
dimensions while the latter enter as a boundary conditions to the evolution equation. Moreover, 
it was established basing on the resummation of the particular class of Feynman diagrams - 
fermion vacuum polarization bubbles — that the leading ambiguity in the perturbative series 
could be only A 2 /Q 2 . However, the latter approach does not take into account the non-abelian 
nature of QCD and, therefore, could miss some important features since it was guessed that it is 
precisely the diagrams with three-gluon vertices [0 which might be important in resolving these 
apparent paradox. In the present paper we make a first step and evaluate the Wilson loop in order 
0(a 2 ). Later extending the result to the case of nonzero gluon mass, A 2 , we will be able to find 
the first non-analytical term in the expansion w.r.t. A which might be a signal for the first power 
correction || up to limitations discussed below. However, the results given here are of interest in 
their own right. 

2. Time-loop technique. Following the approach of Korchemsky and Marchesini we can express 
the soft part of the factorized DY cross section || 

Q 2 DY ^Q2 — ~ = °dy^dy (Q 2 ) WdY (z) (1) 

via the Fourier transform, Wdy(z), of the vacuum average of the Wilson loop 

Wd Y (z) = | ^jv°»WD Y {y), Wdy (y) = ^-<0|TrTRexp (ig f dx,B^ |0>, (2) 

with a path Cdy shown in Fig. 0(a). The symbol V stands for path- while T for (anti-)time- 
ordering to be explained below. Here u = ^(1 — z) is a total soft gluon energy. 
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It is well known that a time-like cross section could not be related to the imaginary part of 
any T-product of currents. Contrary, it is given by the particular absorptive parts of the Feynman 
diagrams. To be able to extract the imaginary part we are interested in we can label in some 
way the field operators in the amplitudes to the right and to the left of the cut. This can be 
suitably done with the help of the Keldysh-like diagram technique It allows to recast 

the program for the calculation of the particular discontinuities of the Feynman diagrams to the 
operator-like language. Consider, for instance, certain F-matrix element which is given by the 
functional integral M. — f VB flj -B, exp (i f dzC). The cross section of the process is 


a = M ] M = JvB^VBM n b\~ ] n B^ +) exp (i J dz £ (+) - i J dz 


( 3 ) 


Here the ’’plus” and ’’minus” superscripts label the fields from the direct and conjugated ampli¬ 
tudes, respectively. One can accept that they are the components of a unique operator B M = 
composed from the time- and anti-time-ordered fields. Now the Green function for 
the B'n field is a 2 x 2-matrix constructed from the usual Feynman propagator, its conjugated 
analogue and its discontinuity via the Cutkosky rules for the lines connecting the direct and final 
amplitude^. For the gluon propagator in the Feynman gauge we have (d = 4 — 2e is a space-time 
dimension) 


/ Ae^(0|T {B;(x) 0 B‘(0)} |0} = S‘ l g„ 


-iD^ ++ \q) -iD^ + \q) \ 
iD(+-\q) iD(—\q) ) ’ 


( 4 ) 


with the following form of the components^ B^ ++ ^(g) = [ZT '*(<?)]* = ( q 2 3 + «0) _1 , D^~ + \q) = 

—27Ti5 + (q 2 ). 

Using these conventions we can easily write the Wilson loop (||) as follows 

W'dyO) = T( 0 |TVT{tM[—oo, #H[ 9 , —oo]tW [—oo, 0]J>W[0, —oo]} |0>, (5) 

where T stands for the time-ordering for the (+) and anti-time one for (—)-fields and the Wilson 
lines are given by the following expressions 

<F P i [y, -oo] = P exp (ig J ^ da^B „(oiPi + y)j , 

r—oo 

-ig / da 2 P 2 fi B ll (-a 2 p 2 + y) 

J 0 

Here the two quark momenta are defined in terms of the light-like vectors, p\ = pi+n*, p 2 = p 2 -n , 
which fix different tangents, i.e. — and + directions, on the light cone and they are normalized 
according to n 2 = n* 2 = 0 and nn* = 1. 

2 By definition a (—) field stands to the left of (+) ones. 

3 Here and below the subscript “+” on the distribution, A = 6,9, means the positivity of the energy flow through 
the cut: A +(q 2 ) = A(q 2 )d(qo). 



<F_ P2 [-oo,t/] = P exp 
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Let us note that the peculiarities of the renormalization of the Wilson lines with a cusp 
, [L5[ or/and with a path lying on the light cone 


3 0 


were studied before. It was 


found that they require an additional renormalization and corresponding anomalous dimensions 
play an important role in perturbative QCD. These are the properties which we will extensively 
exploit presently. 

Evaluating the one-loop graphs (Fig. [!](&) and mirror conjugated) we can easily obtain un¬ 
renormalized expression 

W (1) = C F — T ^~ e K -^ E L e , (7) 


7r 


where L = 


(—jyQlJ> 2e2lE + io) with /a 2 being the MS scale parameter. Subtracting the poles 
1/e we get well-known expression 0 

wSkv) = K,w m = c F | 


V i+ « 2 d 

2 2 


i 


m 


( 8 ) 


3. Two-loop results. Here we perform the evaluation of the Wilson loop in the fourth order of 
perturbation theory. Expanding the generating functional ([|) in the perturbative series we obtain 
a set of Feynman diagrams given in Fig. |Ij. Note that diagrams with virtual dressing of Wilson 
line vanish identically in dimensional regularization since the integration over loop momentum is 
given by a scaleless integral. Therefore, from several number of cuts only those survive which are 
given by the product of tree diagrams (up to exception of graph (h)). Moreover, on top of this 
among the remaining graphs we will evaluate only those (they are displayed in Fig. [[]) proportional 
to maximal non-abelian and fermion components, i.e. CfCa and CfTfNj, which is possible due 
to the non-abelian exponentiation theorem |T7[. According to it 


I^dy = 1 + E ^ (n) = exp E 


w 


(n) 


(9) 


72—1 


\n=l 


so that W^ = |(u/ 1 )) 2 + ft/ 2 ) and W W = if/ 1 ). In what follows we evaluate ft/ 2 ). 

3.1. Vacuum polarization diagrams. The vacuum polarization diagram in Fig. [l](c) reads 


VJ 


( 2 ) 

(c) 


= -C F g 2 S e g 2e p l+ p 2 - / da 2 / d,a ] 

•l—oo J— oo 

x / 7^Wrfe- i9((T2P2 - CTlPl+2/) D ( " ) (g)D (++) (g)2Imn H 


-(q), 


( 10 ) 


where S e = e e ^ E ln4?r ) and the imaginary part of the one-loop polarization operator is 


InTI^g) = a s 


C A 


(5 - 3e) - 2T F N f {l - e) 


T(2-e) 

J T(4 - 2e) 


' ji 2 e lE 


0+(q 2 )(q 2 g l ,u-q^)- (H) 


After Fourier transformation with the help of formula {x±_ = 0) 


d d q iqx 9+(q 2 ) 


(2tt) 


0 q 2 


= 7T 


- 1 - 


d/2^ 


-1-2 n 


T(1 - n)V{d/2 - n) 

[—2(x_ — -i0)(a; + — 70)] rf / 2_ 


( 12 ) 
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Figure 1: Integration path in the Wilson loop in Eq. (|2|) is shown in (a). Non-zero one-loop graph 
(6). Different topologies of nonvanishing Feynman diagrams contributing to the Wilson loop IFdy 
at order 0(a 2 ) with maximal non-abclian factor in (c) — ( h ). To complete the set one has to add 
mirror symmetrical graphs which can be taken into account by appropriate multiplicity factors, 
m a , for graphs displayed here: nib — mc = m d = Wf — 2, m e = 1, m g = m^ = 4. Full blob stands 
for the sum of the vacuum polarization bubbles due to fermions, gluons and ghosts. 


we get 


«!> = Cf (^) [CU(5 - 3e) - iT F N,(\ - £ )] 


r(l-e)r(2-e)r(l-2e) 
16e 3 r(4 - 2e) 


1 + e 


-2e'y E £ / 


2e 


(13) 

The role of diagram (d) is to cancel the second term in the curly brackets in Eq. (|F|) which 
comes from the “gauge dependent” g^qy-piece of the Landau gluon propagator. Multiplying these 
contributions by appropriate multiplicity factors coming from adding of analogous graphs with 
different attachments of gluon lines we obtain 


w ( 2) _ C F (^y [C A (5 - 3e) - AT F N f ( 1 - e)] 


r(l-e)r(2-e)r(l-2e) 


8e 3 T(4 - 2e) 


e 2e ^L 2e 


(14) 


3.2. Box-type diagrams. Among a number of them most vanish due to the light-like character of 
the paths and because planar graphs are proportional to the fermion Casimir operator squared, C F , 
which is omitted due to non-abelian exponentiation theorem. The evaluation of surviving graphs 
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is the most straightforward in the coordinate space where the propagators look like: D^ ++ \x) = 
_i 7r -< i / 2 r((i/2 — 1)[— x 2 +i0] 1-d//2 , and D(~ + \x) = —|7r _d / 2 r(d/2 — 1)[—2(x_ — *0)(x + — i0)] 1-d//2 
(for x± = 0). Fig. |l|(e) gives[] 


w 


,v I r —oo /-u r —oo /*u 

) = -CVCUg 4 <S 2 e// e (pi+P2 -) 2 / den / dcr 2 / do-2 / dcr( 

z JO J— oo J 0 J —oo 


(15) 


X Zd ( +) (<r£p 2 - OlPi + y)dl { +) (cr 2 p 2 “ + J/) = ~C F C A 


r 2 (i 


8e 4 


-e~ 2eiE L 2e . 


Contribution of Fig. [!](/) is w!V, = . Assembling everything together we have 


w 


( 2 ) _ 


= -C F C A 


a < 


7T 


3 r 2 (i 


16 


-e~ 2nE L 2e . 


(16) 


3.3. Non-abelian diagrams. Finally, let us consider the non-abelian diagrams. The typical contri¬ 
bution, for instance for Fig. |l](g) where the only contribution survives when the loop is cut, looks 

like 

d \ 


( 2 ) 

W l) = 


1 /*0 rO no -1 / f) 

-C F C A g A S 2 eia ie p 1+ P2- / da 2 / dai / da( — 

2 d-oo d-oo d-oo \<dcn 


(9(7 i 


x ( -^-e~ ik{a2P2 -< pi+y) D^—\k) [ 1 ^e- i (<- ,T1 ' >qpi D(- + \q)D(- + \k~q). (17) 

d (27T) rf d (27T)“ 


A 


(2vr) d J ( 2? 

The internal integral can be evaluated according to formula^ 

d d q 


/ (£y e "‘ 8 '- D< " +)( « )£ ’ ( " +,(fc -« ) = -j J pry( 4 n 2 - ,i (*: 2 A 2 0 + (*: 2 ) 

x e -^kx+^(kxp-k 2 x 2 ^2 — 


(18) 


- 1 
2 


— % 


{kx) 2 — k 2 x 2 


with x 2 = 0 and the second one with the help of Eq. ([12) making use of the integral representation 
of confluent hypergeometric function ^F\ 0. Thus we get 


w{l] = C F C A ( — 


7r 


r(l-2e)r 2 (l-e) __ 

64e 4 T(2 - 2e) 


2e7 b 


(19) 


4 Note that the colour factor of the box-type diagrams is C F (C F — ^r) but we have kept only maximally 
non-abelian component. 

5 Let us add an interesting side-remark concerning the evaluation of this integral. Actually, there is no need to 
perform its explicit calculation with exponential weight. The dependence on x can be fixed using the properties of 
traceless tensors alone. Expanding the factor e~ lqx in Taylor series and using the fact that the tensor ... q^ n 

is traceless due to presence of S(q 2 ) in the integrand we can parametrize the integral in terms of the only remaining 
vector k according to ({q^q^ ■■•<3y* n )) = {k^k^ 2 (with {...} standing for symmetrization and trace 

subtraction) and perform tensor contraction with the help of the d-dimensional generalization of Nachtmann’s |l7 


convolution {k^k^ ... k lin }x i 




• • • x,, = 


r(w+i)r(A) 

r(n+A) 


( j o 2 \ ^/ .. \ 7 

C n ( with a = | - 1 and C* being the 


Gegenbauer polynomial [[L9|. Using third generating function for the latter |1S|1 we can sum the series back. An 
rr-free function f2 can be easily evaluated then which results into Eq. 
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Now let us observe a fact that the expressions for diagrams we have considered so far with cut 
gluon Green functions coincide with corresponding contributions evaluated with ordinary Feynman 
propagators. This is not accidental but it is a mere consequence of the fact that their collinear 
and cusp singularities — the only ones we are interested in — appear when the gluons propagate 
along the light-cone directions or they are all on short distances. In these cases cut and ordinary 
propagators coincide. 

There are two contributions coming from different cuts of diagram (h), namely, cut-I corre¬ 
sponding to the propagator combination D[ ++ ^ D 2 ++ ^and cut-II with d[ D 2 +) ^3 By 
the reason stated above we have calculated instead the virtual graph which gives 



C F C 


F^A 



r 2 (i-6) 

64e 4 (l + e) 


1 + e - 2 eF 2 


1,1 + e, —2e 
2 + e, 1 - 2e 



-2e7 B ^2e 


( 20 ) 


with F 2 being Appel function [f20| 

Summing Eqs. ( P^ , PD| ) multiplied by a factor of 4 we obtain . 

3-4- Counter-terms. Before the subtraction of the overall divergences we have to handle the sub¬ 
divergences. For this reason we mention that although the diagrams with virtual subgraphs (gluon 
vertex and propagator corrections) vanish for the case at hand, corresponding counter-terms do 
not. Namely, combining the contributions for uncut vacuum blobs and vertex functions results 
into addendum 


w (2) = A) 1^(1) ( 21 ) 

7T 4 e v ' 

where (3 0 = ^T f Nj — ^Ca is the first expansion coefficients of the QCD /3-function (3(g)/g = fjAr 
This completes the set of non-vanishing contributions we have to analyze. 

4. Evolution equation in two-loop approximation. Now with the results derived in the previous sec¬ 
tions we can find the renormalized expression of the two-loop Wilson loop for the DY production^ 
(■u/ 2) = w^ + w 


w< dy( u) = E, e w^ = ^ ju^ In 3 L + w^ hi 2 L + vjf ] In l| , (22) 

with 

(2) l n „ (2) „ r> /67 C(2)\ 5 

w 3 - ~-^C F f3o, w 2 -C f C a I —- —j-—C F T F Nf, 

«.< 2) = C F C A (^ _ I C( 3)) - M. C F T F N,. (23) 


Note the absence of ln 4 L terms [I], |3] due to famous Lee-Nauenberg-Kinoshita theorem [ P| . 
Moreover, this equation is in complete agreement with explicit two-loop calculation of the ordinary 


5 Here we have used an expansion F 2 ^ 


1,1+e, —2e 
2+e,l —2e 


b 1 ) 


i, 11 = (1 - 2C(2)e 2 - 14C(3)e 3 + 0(e 4 )). 
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Feynman graphs by van Neerven et al. p| (cf. Eq. (4.9) there adding charge renormalization 
counter-term) provided we transform their result to the language of moments so that we receive 
an addendum oc £(2) which makes the results coincide. The fact that the coefficient of the leading 
log, w 3 , is proportional to the /3-function suggests that IFdy satisfies a renormalization group 
equation p, [HJ. With the expression (j22|) at hand it is easy to verify that ITdy( 2 /) indeed respects 
the following evolution equation (In fFoYpy) = w[|y(?/) + w|,y(|/)) 


vS- + @(g)-iz \ 11 'Uvfyj = 2r cmp ( g ) in L + r DY (g), 


(24) 


with r, 


cusp 


dp, dg 

(v) C F + (^) 2 Cf (C a (% - - §7>A/) being well-known universal cusp 


anomalous dimension JT3 and a new process-dependent entry 


r DY (g) = (- 

7T 


Cr (^(^->->>) + ( 


C(2) 


I 


(25) 


where ^-dependent part is in agreement with calculation of Ref. [|6j. An unusual feature of 
Eq. ( |2d) ) is that the “anomalous dimension” 7 = 2r cusp In L + T DY depends explicitly on the 
renormalization scale, p. This means the absence of the multiplicative renormalizability of the 
light-like Wilson line. 

Going to moments of the Fourier transformed Wilson loop, Wdy (q^q) = fo dzz N ~ 1 Wny(z), 
results to a mere substitution of yo in the saddle point approximation by y 0 = — z2-^. The solution 


of the evolution equation (|24l) is given by: 
' 11 N 


W'DY 


QN 0 


— VFdy [Q 


x exp 


r 

’q 2 f P 


— (r 


( p N 2 \ 1 N 

cusp ( a s(p)) hi [ "772 1 F^ DY ( a «(d)) 


(26) 


,Q 2 Nl 


with Nq = e~ lE and the upper limit y which sets the boundary for the maximal soft gluon energy. 

The knowledge of Tdy at order 0(a 2 ) allows to predict DY inclusive production in next-to- 
next-to-leading order approximation in large logarithms provided we have three-loop expression for 
rcusp which can be obtained by explicit calculation of the Wilson line with a cusp or be extracted 
as a coefficient in front of 1/[1 — x] + -distribution of the three-loop quark-quark splitting kernel 
once being evaluated in future. 

5. Conclusion and outlook. In present letter we have calculated two-loop approximation of the 
Wilson loop for the DY process. As we have mentioned in the introduction it is straightforward 
to generalize present consideration by taking into account the finite gluon mass propagator and 
look for non-analyticity in A as a trace of power corrections. This could be considered as an 
attempt to resolve the problem with l/Q-power behaviour in DY reaction — a question which 
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has attracted a lot of attention recently BBS- Bnt one has to be careful with this since even 
in the one-loop calculations there is one-to-one correspondence between finite gluon mass and 
renormalon based schemes || only for the sufficiently inclusive quantities while they give different 
output for observables with weighted final state, e.g. inclusive production in e + e”-annihilation 
[Z2f . More clear signal for the leading power correction would be the first renormalon ambiguity in 
the standard resummation of fermion bubble chains. But the technique for handling the diagrams 
with non-abelian vertex still has to be developed. 


We deeply indebted to G.P. Korchemsky for numerous illuminating discussions and to D.I. 
Kazakov for correspondence on a question related to present study. The author was supported by 
the Alexander von Humboldt Foundation and partially by Russian Foundation for Fundamental 
Research, grant N 96-02-17631. 
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